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METHOD AND SYSTEM FOR RECTILINE ARIZING AN IMAGE OF A 
MICROARRAY HAVING A NON-RECTILINEAR FEATURE 

ARRANGEMENT 

5 Embodiments of the present invention are related to data extraction from 

images of microarrays and, in particular, to a general method and system for 
rectilinearizing a double-density microarray having a non-rectilinear, outermost, 
feature-position arrangement. 

1 0 BACKGROUND OF THE INVENTION 

The present invention is related to microarrays. In order to facilitate 
discussion of the present invention, a general background for microarrays and 
examples of their use is provided below. In the following discussion, the terms 
"microarray," "molecular array," and "array" are used interchangeably. The terms 

15 "microarray" and "molecular array" are well known and well understood in the 
scientific community. As discussed below, a microarray is a precisely manufactured 
tool which may be used in research, diagnostic testing, or various other analytical 
techniques. 

Microarray technologies have gained prominence in biological research and in 
20 diagnostics. Currently, microarray techniques are most often used to determine the 
concentrations of particular nucleic-acid polymers in complex sample solutions. 
Microarray-based analytical techniques are not, however, restricted to analysis of 
nucleic acid solutions, but may be employed to analyze complex solutions of any type 
of molecule that can be optically or radiometrically scanned and that can bind with 
25 high specificity to complementary molecules synthesized within, or bound to, discrete 
features on the surface of a microarray. Because microarrays are widely used for 
analysis of nucleic acid samples, the following background information on 
microarrays is introduced in the context of analysis of nucleic acid solutions following 
a brief background of nucleic acid chemistry. 
30 Deoxyribonucleic acid ("DNA") and ribonucleic acid ("RNA") are linear 

polymers, each synthesized from four different types of subunit molecules. Figure 1 
illustrates a short DNA polymer 100, called an oligomer, composed of the following 
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subunits: (1) deoxy-adenosine 102; (2) deoxy-thymidine 104; (3) deoxy-cytosine 106; 
and (4) deoxy-guanosine 108. When phosphorylated, subunits of DNA and RNA 
molecules are called "nucleotides" and are linked together through phosphodiester 
bonds 110-115 to form DNA and RNA polymers. A linear DNA molecule, such as 
5 the oligomer shown in Figure 1, has a 5' end 118 and a 3 ! end 120. A DNA polymer 
can be chemically characterized by writing, in sequence from the 5' end to the 3' end, 
the single letter abbreviations A, T, C, and G for the nucleotide subunits that together 
compose the DNA polymer. For example, the oligomer 100 shown in Figure 1 can be 
chemically represented as "ATCG." 

10 The DNA polymers that contain the organization information for living 

organisms occur in the nuclei of cells in pairs, forming double-stranded DNA helixes. 
One polymer of the pair is laid out in a 5* to 3* direction, and the other polymer of the 
pair is laid out in a 3' to 5' direction, or, in other words, the two strands are anti- 
parallel. The two DNA polymers, or strands, within a double-stranded DNA helix are 

15 bound to each other through attractive forces including hydrophobic interactions 
between stacked purine and pyrimidine bases and hydrogen bonding between purine 
and pyrimidine bases, the attractive forces emphasized by conformational constraints 
of DNA polymers. Because of a number of chemical and topographic constraints, 
double-stranded DNA helices are most stable when deoxy-adenylate subunits of one 

20 strand hydrogen bond to deoxy-thymidylate subunits of the other strand, and deoxy- 
guanylate subunits of one strand hydrogen bond to corresponding deoxy-cytidilate 
subunits of the other strand. Figures 2A-B illustrates the hydrogen bonding between 
the purine and pyrimidine bases of two anti-parallel DNA strands. AT and GC base 
pairs, illustrated in Figures 2A-B, are known as Watson-Crick ("WC") base pairs. 

25 Two DNA strands linked together by hydrogen bonds forms the familiar helix 
structure of a double-stranded DNA helix. Figure 3 illustrates a short section of a 
DNA double helix 300 comprising a first strand 302 and a second, anti-parallel 
strand 304. Although deoxy-guanylate subunits of one strand are generally paired 
with deoxy-cytidilate subunits from the other strand, and deoxy-thymidilate subunits 

30 in one strand are generally paired with deoxy-adenylate subunits from the other 
strand, non-WC base pairings may occur within double-stranded DNA. 
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Double-stranded DNA may be denatured, or converted into single stranded 
DNA, by changing the ionic strength of the solution containing the double-stranded 
DNA or by raising the temperature of the solution. Single-stranded DNA polymers 
may be renatured, or converted back into DNA duplexes, by reversing the denaturing 
5 conditions, for example by lowering the temperature of the solution containing 
complementary single-stranded DNA polymers. During renaturing or hybridization, 
complementary bases of anti-parallel DNA strands form WC base pairs in a 
cooperative fashion, leading to reannealing of the DNA duplex. 

The ability to denature and renature double-stranded DNA has led to the 

10 development of many extremely powerful and discriminating assay technologies for 
identifying the presence of DNA and RNA polymers having particular base sequences 
or containing particular base subsequences within complex mixtures of different 
nucleic acid polymers, other biopolymers, and inorganic and organic chemical 
compounds. Figures 4 - 7 illustrate the principle of the microarray-based 

15 hybridization assay. A microarray (402 in Figure 4) comprises a substrate upon 
which a regular pattern of features is prepared by various manufacturing processes. 
The microarray 402 in Figure 4, and in subsequent Figures 5-7, has a grid-like 2- 
dimensional pattern of square features, such as feature 404 shown in the top, left-hand 
corner of the microarray. Each feature of the microarray contains a large number of 

20 identical oligonucleotides covalently bound to the surface of the feature. These bound 
oligonucleotides are known as probes. In general, chemically distinct probes are 
bound to the different features of an microarray, so that each feature corresponds to a 
particular nucleotide sequence. 

Once a microarray has been prepared, the microarray may be exposed to a 

25 sample solution of target DNA or RNA molecules (410-413 in Figure 4) labeled with 
fluorophores, chemiluminescent compounds, or radioactive atoms 415-418. Labeled 
target DNA or RNA hybridizes through base pairing interactions to the 
complementary probe DNA, synthesized on the surface of the microarray. Figure 5 
shows a number of such target molecules 502-504 hybridized to complementary 

30 probes 505-507, which are in turn bound to the surface of the microarray 402. 
Targets, such as labeled DNA molecules 508 and 509, that do not contain nucleotide 
sequences complementary to any of the probes bound to microarray surface do not 
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hybridize to generate stable duplexes and, as a result, tend to remain in solution. The 
sample solution is then rinsed from the surface of the microarray, washing away any 
unbound-labeled DNA molecules. In other embodiments, unlabeled target sample is 
allowed to hybridize with the microarray first. Typically, such a target sample has 
5 been modified with a chemical moiety that will react with a second chemical moiety 
in subsequent steps. Then, either before or after a wash step, a solution containing the 
second chemical moiety bound to a label is reacted with the target on the microarray. 
After washing, the microarray is ready for scanning. Biotin and avidin represent an 
example of a pair of chemical moieties that can be utilized for such steps. 

10 As shown in Figure 6, the bound labeled DNA molecules are detected via 

optical or radiometric scanning. Optical scanning involves exciting labels of bound 
labeled DNA molecules with electromagnetic radiation of appropriate frequency and 
detecting fluorescent emissions from the labels, or detecting light emitted from 
chemiluminescent labels. When radioisotope labels are employed, radiometric 

15 scanning can be used to detect the signal emitted from the hybridized features. 
Additional types of signals are also possible, including electrical signals generated by 
electrical properties of bound target molecules, magnetic properties of bound target 
molecules, and other such physical properties of bound target molecules that can 
produce a detectable signal. Optical, radiometric, or other types of scanning produce 

20 an analog or digital representation of the microarray as shown in Figure 7, with 
features to which labeled target molecules are hybridized similar to 706 optically or 
digitally differentiated from those features to which no labeled DNA molecules are 
bound. Features displaying positive signals in the analog or digital representation 
indicate the presence of DNA molecules with complementary nucleotide sequences in 

25 the original sample solution. Moreover, the signal intensity produced by a feature is 
generally related to the amount of labeled DNA bound to the feature, in turn related to 
the concentration, in the sample to which the microarray was exposed, of labeled 
DNA complementary to the oligonucleotide within the feature. 

When a microarray is scanned or otherwise analyzed, data may be collected as 

30 a two-dimensional digital image of the microarray, each pixel of which represents the 
intensity of phosphorescent, fluorescent, chemiluminescent, or radioactive emission 
from an area of the microarray corresponding to the pixel. A microarray data set may 
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comprise a two-dimensional image or a list of numerical or alphanumerical pixel 
intensities, or any of many other computer-readable data sets. An initial series of 
steps employed in processing digital microarray images includes constructing a 
regular coordinate system for the digital image of the microarray by which the 
5 features within the digital image of the microarray can be indexed and located. For 
example, when the features are laid out in a periodic, rectilinear pattern, a rectilinear 
coordinate system is commonly constructed so that the positions of the centers of 
features lie as closely as possible to intersections between horizontal and vertical 
gridlines of the rectilinear coordinate system, alternatively, exactly half-way between 

10 a pair of adjacent horizontal and a pair of adjacent vertical grid lines. Then, regions 
of interest ("ROIs") are computed, based on the initially estimated positions of the 
features in the coordinate grid, and centroids for the ROIs are computed in order to 
refine the positions of the features. Once the position of a feature is refined, feature 
pixels can be differentiated from background pixels within the ROI, and the signal 

1 5 corresponding to the feature can then be computed by integrating the intensity over 
the feature pixels. 

Scanning of a microarray by an optical scanning device or radiometric 
scanning device generally produces an image comprising a rectilinear grid of pixels, 
with each pixel having a corresponding signal intensity. These signal intensities are 

20 processed by an microarray-data-processing program that analyzes data scanned from 
an microarray to produce experimental or diagnostic results which are stored in a 
computer-readable medium, transferred to an intercommunicating entity via electronic 
signals, printed in a human-readable format, or otherwise made available for further 
use. Microarray experiments can indicate precise gene-expression responses of 

25 organisms to drugs, other chemical and biological substances, environmental factors, 
and other effects. Microarray experiments can also be used to diagnose disease, for 
gene sequencing, and for analytical chemistry. Processing of microarray data can 
produce detailed chemical and biological analyses, disease diagnoses, and other 
information that can be stored in a computer-readable medium, transferred to an 

30 intercommunicating entity via electronic signals, printed in a human-readable format, 
or otherwise made available for further use. 
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The microarrays illustrated in Figures 4-7 show a regularly spaced two- 
dimensional matrix or microarray features. However, the microarray features may be 
arranged to cover the surface of the microarray at higher densities, as, for example, by 
offsetting the features in adjacent rows to produce a more densely-packed 
5 arrangement of features. Figures 8A-B illustrate two types of arrangements of 
densely-packed, disk-shaped microarray features. Figure 8A shows a densely-packed 
microarray having a rectangular, outermost, feature-position arrangement denoted by 
shaded features 801. Figure 8B shows a densely-packed microarray having a non- 
rectangular, outermost, feature-position arrangement denoted by shaded features 802. 
10 A microarray having a rectangular, outermost, feature-position arrangement is 
referred to as "rectilinear," and a microarray having a non-rectangular, outermost, 
feature-position arrangement is referred to as "non-rectilinear." 

Many currently available programs that extract feature information from 
scanned or otherwise analyzed images of microarrays rely on a rectilinear, outermost, 
1 5 feature-position arrangement in order to generate an initial, or seed, location for each 
feature. Attempts to extract data from microarrays with non-rectilinear, outermost 
feature boundaries by current feature extraction programs often result in errors. 
Therefore, designers, manufacturers, and users of microarrays have recognized a need 
for a method and system to enable extraction of data from microarrays having non- 
20 rectilinear, outermost feature boundaries. 



SUMMARY OF THE INVENTION 

One embodiment of the present invention comprises a method and system for 
rectilinearizing a double-density, non-rectilinear microarray of features within a 

25 scanned or otherwise analyzed image of a microarray. A feature-coordinate grid of 
horizontal and vertical grid lines is superimposed over the microarray image so that 
the center of each feature of the microarray image coincides with a unique 
intersection of a horizontal and vertical grid line. Three corner features are selected 
and indexed. The coordinates of the three selected corner features are used to 

30 determine three feature positions defining three corners of a rectilinear, outermost, 
feature-position arrangement of the non-rectilinear microarray of features. A fourth 



Docket No. 10030772-1 



feature position of the rectilinear, outermost, feature-position arrangement is 
determined from two of the three feature positions. 

BRIEF DESCRIPTION OF THE DRAWINGS 
5 Figure 1 illustrates a short DNA polymer. 

Figures 2A-B illustrate the hydrogen bonding between the purine and 
pyrimidine bases of two anti-parallel DNA strands. 

Figure 3 illustrates a short section of a DNA double helix comprising a first 
strand and a second, anti-parallel strand. 
10 Figure 4 illustrates a grid-like, two-dimensional pattern of square features. 

Figure 5 shows a number of target molecules hybridized to complementary 
probes, which are in turn bound to the surface of the microarray. 

Figure 6 illustrates the bound labeled DNA molecules detected via optical or 
radiometric scanning. 

1 5 Figure 7 illustrates optical, radiometric, or other types of scanning produced 

by an analog or digital representation of the microarray. 

Figures 8A-B illustrate two types of arrangements of densely-packed, disk- 
shaped microarray features. 

Figure 9 is a control-flow diagram of a rectilinearization method that 
20 represents one embodiment of the present invention. 

Figure 10 illustrates the orientation of a double-density, hypothetical 
microarray having a non-rectilinear, outermost feature arrangement. 

Figure 11 shows the two rectilinear lattices of features of the hypothetical 
microarray shown in Figure 10. 
25 Figure 12 shows three example selected corner features of a hypothetical 

microarray. 

Figure 13 illustrates the assignment of coordinates to the features of a scanned 
or otherwise analyzed image of a microarray that represents one embodiment of the 
present invention. 

30 Figure 14 shows the selected corner features F Corner index and the corresponding 

points L x C CornerJndex and L 2 C Corner _ index of a hypothetical microarray. 
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Figure 15 is an illustration of the selected corner features F Comer index9 the 
points L x C Corner Jndex and L 2 C Corner indexi and lattices L/ and L 2 of a hypothetical 
microarray. 

Figure 16 is an illustration of three of the four feature positions defining a 
5 rectilinear, outermost, feature-position arrangement for a hypothetical microarray. 

Figure 17 is an illustration of two orthogonal vectors defining a rectilinear, 
outermost, feature-position arrangement of a hypothetical microarray. 

DETAILED DESCRIPTION OF THE INVENTION 

10 The present invention is directed toward a method and system for 

rectilinearizing a scanned or otherwise analyzed image of a non-rectilinear grid of 
microarray features. Figure 9 is a control-flow diagram of a rectilinearization method 
that represents one embodiment of the present invention. First, in step 901, an image 
of a microarray having a non-rectilinear, outermost, feature-position arrangement is 

15 provided. Next, in step 902, three corner features of the microarray of features are 
selected. Next, in step 903, the coordinates of the three selected corner features are 
used to determine three feature positions defining two orthogonal sides of a 
rectilinear, outermost, feature-position arrangement of the microarray features. 
Finally, in step 904, a fourth feature position is determined from the coordinates of 

20 two of the three feature positions determined in step 903 to complete the rectilinear, 
outermost, feature-position arrangement. 

The following discussion includes two subsections, a first subsection including 
additional information about molecular arrays, and a second subsection describing 
embodiments of the present invention with reference to Figures 10 - 17. 

25 

Additional Information About Microarravs 

An array may include any one-, two- or three-dimensional arrangement of 
addressable regions, or features, each bearing a particular chemical moiety or 
30 moieties, such as biopolymers, associated with that region. Any given array substrate 
may carry one, two, or four or more arrays disposed on a front surface of the 
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substrate. Depending upon the use, any or all of the arrays may be the same or 
different from one another and each may contain multiple spots or features. A typical 
array may contain more than ten, more than one hundred, more than one thousand, 
more ten thousand features, or even more than one hundred thousand features, in an 
5 area of less than 20 cm 2 or even less than 10 cm 2 . For example, square features may 
have widths, or round feature may have diameters, in the range from a 10 (xm to 1.0 
cm. In other embodiments each feature may have a width or diameter in the range of 
1.0 jim to 1.0 mm, usually 5.0 jxm to 500 (im, and more usually 10 jam to 200 |jm. 
Features other than round or square may have area ranges equivalent to that of 

10 circular features with the foregoing diameter ranges. At least some, or all, of the 
features may be of different compositions (for example, when any repeats of each 
feature composition are excluded the remaining features may account for at least 5%, 
10%, or 20% of the total number of features). Inter- feature areas are typically, but not 
necessarily, present. Inter-feature areas generally do not carry probe molecules. Such 

15 inter-feature areas typically are present where the arrays are formed by processes 
involving drop deposition of reagents, but may not be present when, for example, 
photolithographic array fabrication processes are used. When present, interfeature 
areas can be of various sizes and configurations. 

Each array may cover an area of less than 100 cm 2 , or even less than 50 cm 2 , 

20 10 cm 2 or 1 cm 2 . In many embodiments, the substrate carrying the one or more arrays 
will be shaped generally as a rectangular solid having a length of more than 4 mm and 
less than 1 m, usually more than 4 mm and less than 600 mm, more usually less than 
400 mm; a width of more than 4 mm and less than 1 m, usually less than 500 mm and 
more usually less than 400 mm; and a thickness of more than 0.01 mm and less than 

25 5.0 mm, usually more than 0.1 mm and less than 2 mm and more usually more than 
0.2 and less than 1 mm. Other shapes are possible, as well. With arrays that are read 
by detecting fluorescence, the substrate may be of a material that emits low 
fluorescence upon illumination with the excitation light. Additionally in this 
situation, the substrate may be relatively transparent to reduce the absorption of the 

30 incident illuminating laser light and subsequent heating if the focused laser beam 
travels too slowly over a region. For example, a substrate may transmit at least 20%, 
or 50% (or even at least 70%, 90%, or 95%), of the illuminating light incident on the 
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front as may be measured across the entire integrated spectrum of such illuminating 
light or alternatively at 532 nm or 633 nm. 

Arrays can be fabricated using drop deposition from pulsejets of either 
polynucleotide precursor units (such as monomers) in the case of in situ fabrication, 
5 or the previously obtained polynucleotide. Such methods are described in detail in, 
for example, US 6,242,266, US 6,232,072, US 6,180,351, US 6,171,797, US 
6,323,043, U.S. Patent Application Serial No. 09/302,898 filed April 30, 1999 by 
Caren et al., and the references cited therein. Other drop deposition methods can be 
used for fabrication, as previously described herein. Also, instead of drop deposition 

10 methods, photolithographic array fabrication methods may be used. Interfeature areas 
need not be present particularly when the arrays are made by photolithographic 
methods as described in those patents. 

A molecular array is typically exposed to a sample including labeled target 
molecules, or, as mentioned above, to a sample including unlabeled target molecules 

15 followed by exposure to labeled molecules that bind to unlabeled target molecules 
bound to the array, and the array is then read. Reading of the array may be 
accomplished by illuminating the array and reading the location and intensity of 
resulting fluorescence at multiple regions on each feature of the array. For example, a 
scanner may be used for this purpose, which is similar to the AGILENT 

20 MICROARRAY SCANNER manufactured by Agilent Technologies, Palo Alto, CA. 
Other suitable apparatus and methods are described in published U.S. patent 
applications 20030160183A1, 20020160369A1, 20040023224A1, and 
20040021055A, as well as U.S. patent 6,406,849. However, arrays may be read by 
any other method or apparatus than the foregoing, with other reading methods 

25 including other optical techniques, such as detecting chemiluminescent or 
electroluminescent labels, or electrical techniques, for where each feature is provided 
with an electrode to detect hybridization at that feature in a manner disclosed in US 
6,251,685, and elsewhere. 

A result obtained from reading an array, followed by application of a method 

30 of the present invention, may be used in that form or may be further processed to 
generate a result such as that obtained by forming conclusions based on the pattern 
read from the array, such as whether or not a particular target sequence may have 
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been present in the sample, or whether or not a pattern indicates a particular condition 
of an organism from which the sample came. A result of the reading, whether further 
processed or not, may be forwarded, such as by communication, to a remote location 
if desired, and received there for further use, such as for further processing. When 
5 one item is indicated as being remote from another, this is referenced that the two 
items are at least in different buildings, and may be at least one mile, ten miles, or at 
least one hundred miles apart. Communicating information references transmitting 
the data representing that information as electrical signals over a suitable 
communication channel, for example, over a private or public network. Forwarding 

10 an item refers to any means of getting the item from one location to the next, whether 
by physically transporting that item or, in the case of data, physically transporting a 
medium carrying the data or communicating the data. 

As pointed out above, array-based assays can involve other types of 
biopolymers, synthetic polymers, and other types of chemical entities. A biopolymer 

15 is a polymer of one or more types of repeating units. Biopolymers are typically found 
in biological systems and particularly include polysaccharides, peptides, and 
polynucleotides, as well as their analogs such as those compounds composed of, or 
containing, amino acid analogs or non-amino-acid groups, or nucleotide analogs or 
non-nucleotide groups. This includes polynucleotides in which the conventional 

20 backbone has been replaced with a non-naturally occurring or synthetic backbone, and 
nucleic acids, or synthetic or naturally occurring nucleic-acid analogs, in which one or 
more of the conventional bases has been replaced with a natural or synthetic group 
capable of participating in Watson-Crick-type hydrogen bonding interactions. 
Polynucleotides include single or multiple-stranded configurations, where one or 

25 more of the strands may or may not be completely aligned with another. For 
example, a biopolymer includes DNA, RNA, oligonucleotides, and PNA and other 
polynucleotides as described in US 5,948,902 and references cited therein, regardless 
of the source. An oligonucleotide is a nucleotide multimer of about 10 to 100 
nucleotides in length, while a polynucleotide includes a nucleotide multimer having 

30 any number of nucleotides. 

As an example of a non-nucleic-acid-based molecular array, protein antibodies 
may be attached to features of the array that would bind to soluble labeled antigens in 
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a sample solution. Many other types of chemical assays may be facilitated by array 
technologies. For example, polysaccharides, glycoproteins, synthetic copolymers, 
including block copolymers, biopolymer-like polymers with synthetic or derivitized 
monomers or monomer linkages, and many other types of chemical or biochemical 
5 entities may serve as probe and target molecules for array-based analysis. A 
fundamental principle upon which arrays are based is that of specific recognition, by 
probe molecules affixed to the array, of target molecules, whether by sequence- 
mediated binding affinities, binding affinities based on conformational or topological 
properties of probe and target molecules, or binding affinities based on spatial 

10 distribution of electrical charge on the surfaces of target and probe molecules. 

Scanning of a molecular array by an optical scanning device or radiometric 
scanning device generally produces an image comprising a rectilinear grid of pixels, 
with each pixel having a corresponding signal intensity. These signal intensities are 
processed by an array-data-processing program that analyzes data scanned from an 

15 array to produce experimental or diagnostic results which are stored in a computer- 
readable medium, transferred to an intercommunicating entity via electronic signals, 
printed in a human-readable format, or otherwise made available for further use. 
Molecular array experiments can indicate precise gene-expression responses of 
organisms to drugs, other chemical and biological substances, environmental factors, 

20 and other effects. Molecular array experiments can also be used to diagnose disease, 
for gene sequencing, and for analytical chemistry. Processing of molecular-array data 
can produce detailed chemical and biological analyses, disease diagnoses, and other 
information that can be stored in a computer-readable medium, transferred to an 
intercommunicating entity via electronic signals, printed in a human-readable format, 

25 or otherwise made available for further use. 

Embodiments of the Present Invention 

One embodiment of the method of the present invention is described by 
30 applying the method to an example hypothetical microarray having a double-density, 
non-rectilinear, outermost, feature-position arrangement. Figure 10 illustrates the 
orientation of the densely-packed, hypothetical microarray 1002 having a non- 
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rectilinear, outermost, feature-position arrangement. The densely-packed microarray 
shown in Figure 10 may also be referred to as a "double-density" microarray. The 
non-rectilinear, outermost, feature-position arrangement of the hypothetical 
microarray 1002 is denoted by shaded features 1004. The microarray 1002 is oriented 
5 so that corner features 1006 and 1008 of the outermost, feature-position arrangement 
of features 1004 are located in the top, left-hand corner and the bottom, right-hand 
corner of the microarray, respectively. 

Three features of the microarray are selected from the corner features of two 
rectilinear lattices of the hypothetical microarray 1002. Figure 11 shows the two 

10 rectilinear lattices 1101 and 1102, referred to as "Z,y," and "L^, respectively, of the 
hypothetical microarray 1002. The three features are selected from the corner 
features 1 103 - 1 106 of lattice Lj 1101 and the corner features 1 107 - 1 1 10 of lattice 
L2 1102. The corner features are determined by choosing one feature from each of 
three of the four pairs 1111 - 1114 of two corner features. For example, the first 

15 corner feature can be selected from any one of the four pairs 1111 - 1114. The 
second corner feature can be selected from any of the three pairs not already used in 
the selection of the first corner feature. Finally, the third corner feature can be 
selected from the two remaining pairs not already used in the selection of the first and 
second corner features. 

20 After the three corner features have been selected, the selected corner features 

are indexed 1, 2, 3, or 4. In the present example, the indexes are assigned in a 
clockwise manner beginning with assigning the index value 1 to the selected corner 
feature located in the top, left-hand corner the microarray 1002, as denoted by pair 
1111, and ending with assigning the index value 4 to the selected corner feature 

25 located in the bottom, left-hand corner of the microarray, as denoted by pair 1114. 
The selected corner features are referred to as " F Corner index y" where the subscript 

Corner index is assigned the integer value 1, 2, 3, or 4 associated with the corners of 
the lattices Lj and L2. Figure 12 shows three example selected corner features F 4 
1201, Fj 1202, and F 2 1202 of the hypothetical microarray 1002. Note that the 
30 present invention is not limited to selecting the three selected corner features having 
the subscript Corner Jndex equal to 4, 2, 1, as shown in Figure 12. Different possible 
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selected corner features for Corner Jndex equal to 1, 2, 3; or 2, 3, 4; or 4, 3, 1 can also 
be selected. 

In order to rectilinearize a microarray having a non-rectilinear, outermost, 
feature-position arrangement, a feature-coordinate grid identifying the centroid 
5 coordinates of each feature of the microarray is needed. Figure 13 illustrates the 
assignment of coordinates, to the features of an image of a microarray, that represents 
one of many possible embodiments of the present invention. Figure 13 illustrates the 
assignment of coordinates of a feature-coordinate grid to an image of the hypothetical 
microarray 1002. A set of horizontal and vertical grid lines, such as horizontal line 

10 1303 and vertical line 1304, are superimposed on the image of the microarray so that 
the vertical and horizontal grid line intersections coincide as closely as possible with 
the centers of microarray features to give a regularly spaced feature-coordinate grid. 
The distance 1305 between the vertical grid lines is referred to as "offset_x," and the 
distance 1306 between the horizontal grid lines is referred to as "offset_y." The 

15 feature-coordinate grid establishes a two-dimensional system for specifying the 
location of each microarray feature. Thus, for example, the location of feature 1302 
can be specified by the grid-point indices (2,0). 

Next, the x and y coordinates of each of the three selected corner features 
F comer jndex are usec * to determine the coordinates of the two points L x C Comer index and 

20 L 2 C Corner index for each of the selected corner features according to the following 
expressions: 

Lfi Corner Jndex — X ~ & Corner Jndex — X ~ °ff Set _ X 

(1) 

L \ C Corner Jndex _ V = F Comer Jndex _ ^ ~ °ff Set - J> 

(2) 

L 2 C Corner Jndex & Corner Jndex — X ~ °ff Set _ X 

(3) 

L 2 C Corner Jndex - ^ = F Corner Jndex - > " °ff Set - J 7 

(4) 



25 



30 
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where LiCcomerjndex_x and LiCcomerjndex_y are the x and y coordinates of the point 
L \ C comer_index> respectively, L 2 C 'comer jndexjc and L 2 C Comer judex are the x and y 
coordinates of the point L 2 C CortterJndex , respectively, and F Co rnerjndex_x and 
5 Fcomer_index_y are the x and y coordinates of the selected corner feature F Corner index . 

Figure 14 shows the selected corner features F Corner index and the corresponding 
points L x C Corner index and L 2 C CornerJndex determined according to equations (1) through 
(4) for the hypothetical microarray 1002. The centroid coordinates of corner feature 
Fj 1202 are used to determine the points L l C l 1401 and L 2 C/ 1402, as indicated by 
10 arrows 1403 and 1404, respectively. The centroid coordinates of feature F 2 1203 are 
used to determine the points L X C 2 1405 and L 2 C 2 1406, as indicated by arrows 1407 
and 1408, respectively. Finally, the centroid coordinates of corner feature F 4 1201 are 
used to determine the points L X C A 1409 and L 2 C 4 1410, as indicated by edges 1411 
and 1412, respectively. 

15 Next, the selected corner features F Corner index are determined to be either a 

member of lattice Lj or L 2 according the following four conditions: 

(1) For selected corner feature F/, if each of the points LjCj and L 2 Ci coincide 
with features of the microarray, then selected corner feature Fj is a corner feature of 
lattice L 2 . If LjCj does not coincide with a feature of the microarray, and L 2 Cj 

20 coincides with a feature of the microarray, then selected corner feature Fj is a corner 
feature of lattice Lj . 

(2) For selected corner feature F 2 , if the point L/C 2 coincides with a feature of 
the microarray, and L 2 C 2 does not coincides with any features of the microarray, then 
selected corner feature F 2 is a corner feature of lattice L 2 . If L t C 2 does not coincide 

25 with a feature of the microarray, and L 2 C 2 coincides with a feature of the microarray, 
then selected corner feature F 2 is a corner feature of lattice Lj. 

(3) For selected corner feature Fj, if each of the points Z7C3 and L 2 C$ coincide 
with features of the microarray, then selected corner feature F 3 is a corner feature of 
lattice L/. If L1C3 coincides with a feature of the microarray, and L 2 Cs does not 
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coincides with a feature of the microarray, then selected corner feature F3 is a corner 
feature of lattice L2. 

(4) For selected corner feature F4, if the point L1C4 coincides with a feature of 
the microarray grid, and L2C4 does not coincides with any features of the microarray, 
5 then selected corner feature F4 is a corner feature of lattice L2. If L1C4 does not 
coincide with a feature of the microarray, and L2C4 coincides with a feature of the 
microarray, then selected corner feature F 4 is a corner feature of lattice L/. 

Figure 15 is an illustration of the hypothetical selected corner features 

Fcomerjndex* Ae P° intS ^Corner Jndex *"I L 2 C Corner Jndex , and lattices Lj and L 2 Of the 

10 hypothetical microarray 1002. The points LjCj 1401 and L2C1 1402 coincide with 
features 1501 and 1502, therefore, according to condition (1), selected corner feature 
F] 1202 is a corner feature of lattice L 2 1102. The point L1C2 1405 coincides with 
feature 1403, and the point L2C2 1406 does not coincide with any features of the 
microarray, therefore, according to condition (2), the selected corner feature F 2 1203 

15 is a corner feature of lattice L 2 1 102. Finally, the point LjC 4 1409 coincides with the 
feature 1504, and the point L2C4 1420 does not coincide with any feature of the 
microarray, therefore, according to condition (4), the selected corner feature F 4 1201 
is a corner feature of lattice L2 1 102. 

Next, the coordinates of three of the four feature positions defining the four 

20 corner features of the rectilinear, outermost, feature-position arrangement of the 
microarray of features are determined. The feature positions defining the four corner 
features of a derived, outermost, rectilinear arrangement are referred to as 
" p comer jndex" where the values of the subscript Corner index is assigned the values 1, 
2, 3, and 4 that coincide with corner indices as describe above in relation to Figure 1 1 . 

25 The x and y coordinates of three of the four feature positions P Corner index > referred to as 

F Corner index. 

x and P comer jndex_y, respectively, are determined according to following 
eight conditions: 

(1) If selected corner feature Fj is a corner feature of lattice Zy, then 



30 



P/_x = Fj_x 9 and 
Pijy = F/_y. 
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(2) If selected corner feature Fj is a corner feature of lattice L? 9 then 

Pj_x = LjCj x, and 

(3) If selected corner feature F2 is a corner feature of lattice Lj> then 

P2_x = LjC2_x, and 
P 2 _y = F 2 _y. 

(4) If selected corner feature F2 is a corner feature of lattice L 2 > then 

P2_x ~ F2_x, and 
P 2 ^ = LjC 2 jy 

(5) If selected corner feature F3 is a corner feature of lattice Z,/, then 

= Fj_;c, and 
P 3 ^y = F 3 _y. 

(6) If selected corner feature F3 is a corner feature of lattice L 2 , then 

P 3 _x = L 2 C 3 _x 9 and 

(7) If selected corner feature F4 is a corner feature of lattice Lj 9 then 

P*_.x = F4_x 9 and 
P 4 _y = L 2 C4_y. 

(8) If selected corner feature is a corner feature of lattice L 2 > then 
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P4JC = L/C4_x, and 
P4 _y = F 4 _y. 

5 Figure 16 is an illustration of three of the four feature positions defining the 

rectilinear, outermost, feature-position arrangement for the hypothetical microarray 
1002. Because selected corner feature Fj 1202 is a corner feature of lattice L2 1 102, 
then, according to the condition (2), the feature position P/ 1601 is assigned the x and 
y coordinate values of the point L/C/, respectively. Because selected corner feature 

10 F2 1203 is a corner feature of lattice L2 1 102, then, according to the condition (4), the 
feature position P2 1 602 is assigned the x coordinate of selected corner feature F2 and 
the y coordinate of point L/C2. Finally, because selected corner feature F3 1201 is a 
corner feature of lattice L2 1 102, then, according to condition (8), the feature position 
P 2 1604 is assigned the x coordinate of point L1C4 and the y coordinate of selected 

1 5 corner feature F 4 . 

The three feature positions P/ 1601, P2 1602, and P 4 1604 define two vectors 
v/2 and v/4 given by: 

v 12 =P 2 -^=<7,0) 

20 v 14 =/>-/> =<0,5) 

where the brackets "( )" represent vector coordinates. Figure 17 is an illustration of 

the two vectors v/2 1701 and vj 4 1702 defining a rectilinear, outermost, feature- 
position arrangement of the hypothetical microarray 1002. The vectors v l2 1701 and 
25 v/ 4 1702 are orthogonal vectors in the feature-coordinate grid as verified by the scalar 
product given by: 

v l2 -v 14 =(7,0)(0,5> = 0 

30 Because a vector is a directed line segment having an initial point and an ending 
point, the line segments associated with the two vectors v I2 1701 and v i4 1702 form 
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two sides of the derived, rectilinear, outermost, feature-position arrangement of the 
microarray. 

The fourth feature position completing the derived, rectilinear, outermost, 
feature-position arrangement of the microarray is determined according to the 
5 following four conditions: 

(1) If the feature positions Pj, P 2i and P 3 are known, then 



P 4 _x = P } _x 
P4_y = Ps_y 

10 

(2) If the feature positions P 2 , P 3 , and P 4 are known, then 



Pijc = P 4 _x 
Pijy = P 2 jy 

15 

(3) If the feature positions Pj, P 3 , and P 4 are known, then 



P 2 _x = P 3 _x 
P2jy = Pi_y 

20 

(4) If the feature positions Pj, P 2 , and P 4 are known, then 



25 



30 



P 3 _x = P 2 _x 

Ps_y = P4jy 

In Figure 17, the fourth feature position P 3 1703 of the hypothetical 
microarray 1002 is determined according to condition (4). The three feature positions 
P 2 1602, P 3 , 1703, and P 4 1604 define the orthogonal vectors v 23 1704 and v 43 1705 
given by: 

V 2 3=/ > 3-^=(0,5) 
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v 43 =^3-^4 =(7,0) 
where v 12 *v 43 =0 

The vector v 2 3 1704 is orthogonal to the vector v I2 1701, and the vector v 43 1705 is 
5 orthogonal to the vector v I4 1702. The line segments associated with the vectors v 23 
1704 and v 43 1705 provide the two remaining sides of a rectilinear, outermost, feature- 
position arrangement of the microarray 1002. Finally, the two feature positions P 2 
1602 and P 4 1604 are added as imaginary features to the non-rectilinear, outermost, 
feature-position arrangement to give a rectilinear, outermost, feature-position 

10 arrangement of the hypothetical microarray 1002. 

The C-like pseudo code implementation showing an embodiment of the 
present invention is provided below. Note that the pseudo code implementation is not 
intended to describe a complete rectilinearization program for the microarray feature 
data, but to provide sufficient detail to illustrate one possible embodiment of the 

15 rectilinearization methodology as the embodiment might occur within a microarray 
feature extraction program or in microarray feature extraction and data processing 
equipment. The rectilinearization program utilizes a feature-coordinate grid similar to 
that described above in relation to Figure 13 to produce a derived, outermost, 
rectilinear feature-position arrangement of the microarray feature data positions that 

20 are stored in a computer-readable medium, transferred to an intercommunicating 
entity via electronic signals, printed in human-readable format, or otherwise made 
available for further use. 

First, the pseudocode implementation includes several constants: 

25 1 const int offset_x; 

2 const int offset_y ; 

3 const int M ; 

4 const int N ; 

30 The constants "offset_x" and "offset_y" describe the spacing between the vertical and 
horizontal grid lines, respectively, as described above in relation to Figure 13. The 
constants "M" and "N" describes the number of horizontal grid lines, and vertical grid 
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15 



30 



35 



40 
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lines. The specific values for the constants depend on the particular microarray image 
being analyzed. 

Next, the pseudocode implementation includes the structure "xy_coord" 
provided below: 

1 struct xy_coord 

2 { 

3 int x; 

4 inty; 

5 }; 

The variables "x" and "y" are the rectangular Cartesian x and y coordinates used to 
describe coordinate locations in the feature-coordinate grid as described above in 
relation to Figure 13. 

Next, the pseudocode implementation contains the two function prototypes: 



1 void corner_yertices(xy_coord *f,xy_coord *p); 

2 void fourth_point(xy_coord *p,int *ci); 



The functions "corner_vertices" are "fourth_point" are identified as void functions. 
20 The function "corner_vertices" expects the structure pointers "f* and "p" as 

arguments. The function "fourth_point" expects the structure pointer "p" and the 

pointer "ci" as arguments. 

The function "corner_vertices" contains the code relevant to one embodiment 

of the present invention. An implementation of the function "corner vertices" is 
25 provided below: 



1 


void corner_vertices(xy_coord *f,xy_coord *p) 


2 


{ 




3 




int latt1_corn_x,latt1_corn_y; 


4 




int latt2_corn_x,latt2_com_y; 


5 




int corner index; 


6 




intc_i[4]={0,0,0,0}; 


7 




int *ci=c_i; 


8 




int existl ,exist2; 


9 

10 




for (int i=0;i<=2;i++) 


11 




{ 


12 




if(0 <= f->x && f->x <= 2) 


13 




if(0 <= f->y && f->y <= 2) 


14 




{ 


15 




corner index = 1 ; 


16 




cj[0] = 1; 
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17 } 

1 8 if(N-2 <= f->x && f->x <= N) 

1 9 if(0 <= f->y && f->y <= 2) 

20 { 

5 21 cornerjndex = 2; 

22 c_i[1] = 1; 

23 } 

24 if(N-2 <= f->x && f->x <= N) 

25 if(M-2 <= f->y && f->y <= M) 
10 26 { 

27 cornerjndex = 3; 

28 c_i[2] = 1; 

29 } 

30 if(0 <= f->x && f->x <= 2) 

1 5 31 if(M-2 <= f->y && f->y <= M) 

32 { 

33 cornerjndex = 4; 

34 cj[3] = 1; 

35 } 
20 36 

37 latt1_corn_x = f->x - offset_x; 

38 latt1_corn_y = f->y - offset_y; 

39 latt2_corn_x = f->x + offset_x; 

40 latt2_corn_y = f->y + offset_y; 
25 41 

42 exist1=0; 

43 exist2=0; 

44 if (latt1_corn_x >= 0 && latt1_corn_x <= N) 

45 { 

30 46 if (latt1_corn_y >= 0 S.& latt1_corn_y <= M) 

47 existl =1; 

48 } 

49 if (latt2_corn_x >= 0 && latt2_corn_x <= N) 

50 { 

35 51 if (latt2_corn_y >= 0 && latt2_corn_y <= M) 

52 exist2=1; 

53 } 
54 

55 if (cornerjndex -- 1 ) 

40 56 { 

57 if (existl == 0 && exist2 == 1 ) 

58 { 

59 p->x = f->x; 

60 p->y = f->y; 
45 61 } 

62 , if (existl == 1 && exist2 == 1) 

63 { 

64 p->x = latt1_corn_x; 

65 p->y = latt1_corn_y; 
50 66 } 

67 } 

68 if (cornerjndex == 2) 

69 { 
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70 p = P + 1; 

71 if (existl == 0 && exist2 == 1 ) 

72 { 

73 p->x = latt2_corn_x; 
5 74 p->y = f->y; 

75 } 

76 if (existl == 1 && exist2 == 0) 

77 { 

78 p->x = f->x; 

10 79 p->y = latt1_corn_y; 

80 } 

81 P = P-1; 

82 } 

83 if (corner_index -- 3) 
15 84 { 

85 p = p + 2; 

86 if (existl == 1 && exist2 == 0) 

87 { 

88 p->x = f->x; 
20 89 p->y = f->y; 

90 } 

91 if (existl == 1 && exist2 == 1) 

92 { 

93 p->x - latt2_corn_x; 
25 94 p->y = latt2_corn_y; 

95 } 

96 p = p - 2; 

97 } 

98 if (cornerjndex ™ 4) 
30 99 { 

100 p = p + 3; 

101 if (existl == 0 && exist2 == 1) 

102 { 

103 p->x = f->x; 

35 104 p->y = latt2_corn_y; 

105 } 

106 if (existl == 1 && exist2 == 0) 

107 { 

1 08 p->x = latt1_corn_x; 

40 109 p->y = f_>y; 

110 } 

111 P = p-3; 

112 } 

113 f++; 
45 114 } 

1 1 5 fourth_point(p,ci); 

116 } 



The function "corner_vertices" takes the following arguments: (1) "f," a pointer to the 
50 structure elements Fcomerjndexjt and F Corner jndexjy of the array holding the three 
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selected corner features Fc 0 merjndex\ and (2) "p," a pointer to the structure elements 

Pcorner_index_X and P Corner Jndex_y of the array holding the feature positions Pcornerjndex- 

On lines 3 - 8, a number of local variables are declared. These local variables include: 
(1) "lattl_corn_x," "lattl_corn_y," "latt2_corn_x," and "latt2_corn_y," the x and y 
5 coordinates of the points L,C Corner _ index and L 2 C CornerJndex , respectively, as defined 

above in equations (1) - (4); (2) "corner_index," the index of the corner feature of 
lattices Lj and Lj as described above in relation to the hypothetical example shown in 
Figure 11; (3) "c_i," a four-element array, where the element of "c_i" having the 
value "0" identifies the corner point to be determined in the function "fourth_point;" 
10 (4) "ci " a pointer to the elements of array "c_i;" and (5) "existl" and "exist2," 
variables used to determine whether or not the two points L x C Corner index and 

L iC comer jndex coincide with features of the microarray. In the ^br-loop of line 10, lines 
1 1 - 1 14 are repeated for each of the three selected corner features to determine the x 
and y coordinates of three of the feature positions Pcornerjndex- On lines 12 and 13, if 

15 the selected corner feature is located in the top, left-hand corner of either lattice L/ or 
L2, as described above in relation to Figure 11, then on line 15, "corner_index" is 
assigned the value "1," and on line 16, the first element of "c_i" is assigned the value 
"1." On lines 12 and 13, if the selected corner feature is not located in the top, left- 
hand corner of either lattice Lj or £2, then the program proceeds to line 18. On lines 

20 18 - 35, an analogous procedure to that of lines 12 - 17 is carried out to determine the 
value of "cornerindex" and the elements of the array "c_i." On lines 37 - 40, 
"lattl_corn_x," "lattl_corn_y," "latt2_corn_x," and "latt2_corn_y" are calculated, 
respectively. On lines 42 and 43, "existl" and "exist2" are both assigned the value 
"0." On lines 44 and 46, if the point LjCcomerjndex coincides with a feature of the 

25 microarray, then on line 47, "existl" is assigned the value "1." On lines 49 and 51, if 
the point L2Ccomerjndex coincides with a feature of the microarray, then on line 52, 
"exist2" is assigned the value "1." In the nested /^statements of lines 55 - 112, the x 
and y coordinates of three of the four derived, outermost feature positions Pj, P 2 , P3, 
and P 4 are determined according to the conditions 1-8, as described above in relation 

30 to Figure 16. On line 55, if "cornerjndex" equals "1," then lines 57 and 62 are 
executed. On lines 59 and 60, Pj_x and P } _y are assigned the coordinate values F t jc 
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and F/_v, respectively, if the value of "existl" is "0," and the value of "exist2" is "1" 
on line 57. On lines 64 and 65, P/jk and Pi_y are assigned the coordinate values 
LjCj_x and Z,/C/_y, respectively, if the value of "existl" is "1," and the value of 
"exist2" is "1" on line 62. On lines 68 - 82, an analogous procedure to that of lines 
5 55-67 is carried out to determine the structure elements P2__x and P2_y- On lines 83 
- 97, an analogous procedure to that of lines 55 - 67 is carried out to determine the 
structure elements P 3 _x and On lines 98 - 1 12, an analogous procedure to that 

of lines 55 - 67 is carried out to determine the structure elements P 4 _x and P 4 _y. On 
line 1 13, "f ' is incremented. On line 115, the function "fourth_point" is called. 
10 The function "fourth_point" contains the code relevant to one embodiment of 

the present invention. An implementation of the function "fourth_point" is provided 
below: 



I void fourth_point(xy_coord *p,int *ci) 
15 2 { 

3 inttemp_x; 

4 inttemp_y; 

5 if (*ci == 0) 

6 { 

20 7 p = p + 3; 

8 temp_x = p->x; 

9 P = P-2; 

10 temp_y = p->y; 

II p = p-1; 

25 12 p->x = temp_x; 

13 p->y = temp_y; 

14 return; 

15 } 

16 ci++; 

30 17 jf (*ci == 0) 

18 { 

19 p = p + 2; 

20 temp_x = p->x; 

21 P = p-2; 

35 22 temp_y = p->y; 

23 P = P+1; 

24 p->x = temp_x; 

25 p->y = temp__y; 

26 return; 
40 27 } 

28 ci++; 

29 if (*ci == 0) 

30 { 

31 P = P+1; 
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32 temp_x = p->x; 

33 p = p + 2; 

34 temp_y = p->y; 

35 P = P-1; 

5 36 p->x = temp_x; 

37 p->y = temp_y; 

38 return; 

39 } 

40 ci++; 

10 41 if (*ci == 0) 

42 { 

43 temp_x = p->x; 

44 p = p + 2; 

45 temp__y = p->y; 
15 46 P = P+1; 

47 p->x = temp_x; 

48 p->y = temp_y; 

49 return; 

50 } 



20 51 } 

The function "fourth_point" takes the following arguments: (1) "p," a pointer to the 
structure elements PcornerjndexJ* and P Corner jndex_y of the array holding the four, 
outermost feature positions Pj, P 2 , P3, and P 4 defining a rectilinear, outermost feature- 

25 position arrangement of the microarray; and (2) "ci," a pointer to the elements of the 
array "c_i." On lines 3 and 4, two local variables are declared. These local variables 
include: "temp_x" and "temp_y," variables that hold the x and>> coordinate values. If 
the fourth, outermost feature position to be determined is Py, then lines 5 - 15 are 
executed to determine Pijc and P/_y according to condition 1, as described above in 

30 relation to Figure 17. On line 5, if the first element of "c_i" is "0," then lines 7-14 
are executed. On line 7, "p" points to the structure element holding the values P 4 _x 
and P 4 _y. On line 8, "temp_x" is assigned the value of P 4 _x. On line 9, "p" points to 
the structure element holding the values P 2 _x and P 2 _y. On line 10, "temp_y" is 
assigned the value of P 2 _y. On line 11, "p" points to the structure element holding the 

35 values Pj x and Pi_y. On lines 12 and 13, Pj_x and P } _y are assigned the values 
stored in "temp__x" and "temp_y," respectively. On line 14, the program returns to 
the calling function "corner_vertices." On line 5, if the first element of array "c_i" is 
"1," then the program skips lines 6-15. On line 16, "ci" is incremented. On lines 17 
- 27, an analogous procedure to that of lines 5 - 16 is carried out to determine the 

40 structure elements P 2 _x and P 2 _y, if the fourth, outermost feature position to be 
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determined is P2. On lines 29 - 39, an analogous procedure to that of lines 5 - 16 is 
carried out to determine the structure elements Ps_x and P3_y, if the fourth, outermost 
feature position to be determined is P3. On lines 41 - 50, an analogous procedure to 
that of lines 5 - 16 is carried out to determine the structure elements P4JC and if 
5 the fourth, outermost feature position to be determined is P4. 

Although the present invention has been described in terms of a particular 
embodiment, it is not intended that the invention be limited to this embodiment. 
Modifications within the spirit of the invention will be apparent to those skilled in the 
art. For example, an almost limitless number of different implementations of the 

10 many possible embodiments of the method of the present invention can be written in 
any of many different programming languages, embodied in firmware, embodied in 
hardware circuitry, or embodied in a combination of one or more of the firmware, 
hardware, or software, for inclusion in microarray data processing equipment 
employing a computational processing engine to execute software or firmware 

15 instructions encoding techniques of the present invention or including logic circuits 
that embody both a processing engine and instructions. In various embodiments of 
the present invention, the process of selecting corner features can be performed 
automatically by a computer program, or the corner features can be selected manually. 
In various embodiments, the selected corner features can be ordered in various ways 

20 such as in a clockwise or counter-clockwise manner beginning with any one of the 
four corners of the microarray. In various embodiments, different orderings of the 
corner features may be employed. In various embodiments, different possible 
selected comer features may be considered. 

The foregoing description, for purposes of explanation, used specific 

25 nomenclature to provide a thorough understanding of the invention. However, it will 
be apparent to one skilled in the art that the specific details are not required in order to 
practice the invention. The foregoing description of specific embodiments of the 
present invention are presented for purposes of illustration and description. They are 
not intended to be exhaustive or to limit the invention to the precise forms disclosed. 

30 Obviously many modifications and variations are possible in view of the above 
teachings. The embodiments are shown and described in order to best explain the of 
the invention and its practical applications, to thereby enable others skilled in the art 



f 
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to best utilize the invention and various embodiments with various modifications as 
are suited to the particular use contemplated. It is intended that the scope of the 
invention be defined by the following claims and their equivalents: 



